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Collinear splitting, parton evolution and the strange-quark asymmetry of 
the nucleon in NNLO QCD* 
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We consider the collinear limit of QCD amplitudes at one-loop order, and their factorization properties directly 
in colour space. These results apply to the multiple collinear limit of an arbitrary number of QCD partons, and 
are a basic ingredient in many higher-order computations. In particular, we discuss the triple collinear limit and 
its relation to flavour asymmetries in the QCD evolution of parton densities at three loops. As a phenomenological 
consequence of this new effect, and of the fact that the nucleon has non- vanishing quark valence densities, we 
study the perturbative generation of a strange-antistrange asymmetry s{x) — s{x) in the nucleon's sea. 



1. INTRODUCTION 

The high precision of experiments at past, 
present and future particle colliders (LEP, 
HERA, Tevatron, LHC, e+e~ linear colliders) de- 
mands a corresponding precision in theoretical 
predictions. As for perturbative QCD predic- 
tions, this means calculations beyond the next- 
to- leading order (NLO) in the strong coupling as- 
Recent years have witnessed much progress in this 
field. In particular, a great deal of work has been 
devoted to study the properties of QCD scatter- 
ing amplitudes in the infrared (soft and collinear) 
region [T]-[n|. 

The understanding of the infrared singular be- 
haviour of QCD amplitudes is a prerequisite for 
the evaluation of infrared-finite cross sections 
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(and, more generally, infrared- and coUinear-safe 
QCD observables) at higher orders in pertur- 
bation theory. Moreover, the information on 
the infrared properties of the amplitudes can be 
exploited to compute large (logarithmically en- 
hanced) perturbative terms and to resum them 
to all perturbative orders ^Sl- Those investiga- 
tions are also valuable for improving the physics 
content of Monte Carlo event generators (see e.g. 
Ref. 153). In addition, these studies prove to be 
useful even beyond the strict QCD context, and 
can provide hints on the structure of highly sym- 
metric gauge theories at infinite orders in the per- 
turbative expansion (e.g. N=4 super- Yang-Mills, 
see Ref. [Hj). 

Another important application is the calcula- 
tion of the Altarelli-Parisi (AP) kernels, that 
control the scale evolution of parton densities 
and fragmentation functions. The calculation of 
the next-to-next-to- leading order (NNLO) kernels 
has been completed very recently JH|- Collinear 
factorization at the amplitude level (see Sect. [SJ 
can be used jl7| as an alternative and indepen- 
dent method to perform that calculation. To this 
purpose, the one-loop triple collinear splitting |llj 



(see Sect. ^ is one of the necessary ingredients. 
Two other ingredients are the tree-level quadru- 
ple collinear splitting [5] and the two-loop double 
coUinear splitting |12) . 

Besides increasing the quantitative precision 
of the theoretical calculations, the evaluation 
of higher-order contributions can reveal qualita- 
tively new quantum effects. An interesting ex- 
ample is the perturbative generation of charge 
asymmetries in the nucleon's sea TH" (see Sect.0J, 
which arises from the NNLO evolution of parton 
densities. 

2. COLLINEAR FACTORIZATION IN 
COLOUR SPACE 

We consider a generic scattering process in- 
volving final-state QCD partons (massless quarks 
and gluons) of flavour ai,a2,... and momenta 
pi^P2, ■ ■ ., which is described by the matrix ele- 
ment A^ai,a2,...(Pi7P2, ■ ■ •); the external legs are 
on shell (p| = 0) and have physical spin polariza- 
tions. Up to one-loop order, one can write 



M = (gs) 



2n 



0(«s) 



(1) 



where the overall power q is integer. The one-loop 
amplitude M^^' contains ultraviolet and infrared 
singularities that are regularized by using dimen- 
sional regularization in d — 4 — 2e space-time 
dimensions (/i is the dimensional-regularization 
scale). 

The multiple collinear limit is approached when 
the momenta pi , . . . , Pm of m partons become 
parallel. In this limit all the particle subenergies 
Sij = {pi -\- PjY, with i,j = 1, . . . , TO, are of the 
same order and vanish simultaneously, and the 
matrix element ^A{pl, . . . ,pm,Pm+i, ■ ■ ■) becomes 
singular. At the tree level, the dominant sin- 
gular behaviour is M^^\pi,. . . ,Pm,Pm+i, ■ ■ ■) ^ 
{1/y/s)™^^, where s generically denotes a two- 
particle subenergy Sij , or a three-particle suben- 
ergy Sijk, and so forth. At one-loop order, 
this singular behaviour is simply modified by 
scaling violation, M^^\pi, . . . ,Pm,Pm+i, ■ ■ ■) ~ 
(f/-/s)™~^(s/^^)^'^. The dominant singular 
behaviour is captured by universal (process- 
independent) factorization formulae, that are 
usually presented upon decomposition in colour 
subamplitudes |n|~|ni- Collinear factorization is 






Figure f . Factorization of tree-level amplitudes 
in the multiple collinear limit. 
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Figure 2. Factorization of one-loop amplitudes in 
the multiple collinear limit. 

nonetheless valid directly in colour space |11| . 

The colour-space factorization formulae for the 
multiple collinear limit of the tree-level and one- 
loop amplitudes M^^'' and Al'^^-' are: 






,Pm,Pm+l,---)> 
(0) 



(2) 



SpZ..a,SPl, ■ ■ ■ ,Pm) |A^i':L + .,...(^'P™+l' ••■)>, 

l-^S...,a„,a„+i,...bl, • • • ,Prn,Pm+l ■■■))- (3) 

Spi\\..aJPl, ■ ■ ■ ,Pm) \M<;^L^^..iP,Pm+l, ...)> + 
Sp^°l.aJPU ■ ■ ■ ,Pm) |A^i!L + ....(^'P™+l. •••)>• 

These factorization formulae are valid in any 
number d = 4 — 2e of space-time dimensions. 
The only approximation involved on the right- 
hand sides amounts to neglecting terms that are 
less singular in the multiple collinear limit. A 
graphical representation of the factorization for- 
mulae is shown in Figs. Q] and [5] Equations (0) 
and ^ relate the original matrix element (on the 
left-hand side) with m + k partons (where k is ar- 
bitrary) to a matrix element (on the right-hand 
side) with I + k partons. The latter is obtained 



from the former by replacing the m coUinear par- 
tons with a single parent parton, whose momen- 
tum P [P^ — 0) defines the coUinear direction 
and whose flavour a is determined by flavour con- 
servation in the splitting process a — > ai -|- . . . + 
am- The derivation of Eq. |5J in colour space is 
quite straightforward 4 . Its one-loop extension, 
Eq. 0, follows, in particular, from colour coher- 
ence of QCD radiation |19j . 

The process dependence of the factorization 
formulae is entirely embodied in the matrix el- 
ements. The tree-level and one-loop factors 
^Pai---a^ and 5p|jj...a„i which encode the singu- 
lar behaviour in the multiple coUinear limit, are 
universal (process-independent). They depend 
on the momenta and quantum numbers (flavour, 
spin, colour) of the m partons that arise from the 
coUinear splitting. The splitting matrix Sp^ ...^ 
is a matrix in colour+spin space, acting onto the 
colour and spin indices of the m coUinear partons 
on the left and onto the colour and spin indices 
of the parent parton on the right. 

The square of the splitting matrix Sp^^...^^, 
summed over final-state colours and spins and 
averaged over colours and spins of the parent 
parton, defines the m-parton splitting function 
(^ai-a„), which is a generalization of the cus- 
tomary (i.e. with m = 2) AP splitting func- 
tion |2ni- The normalization of the tree-level 
(-Pai--a,„) and one-loop (Pii...a„) splitting func- 
tions is fixed by 
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The one- loop amplitude A^^^^ and, hence, 
Sp^^' have ultraviolet and infrared divergences 
that show up as e-poles in dimensional regular- 
ization. The one-loop splitting matrix can be de- 
composed as 



5pi; 



(1) 



Spi\ 



(l)di- 



Spin 



(5) 



,(l)d 



where S^pii-Xn contains all the e-poles and 
Spi\\^'^;^ is finite when e ^ 0. In Ref. [TT| 



we have presented the explicit expression of 
SPa ■■■a'' for an arbitrary number m of final-state 
coUinear partons. 

3. ONE-LOOP TRIPLE COLLINEAR 
SPLITTING 

The one-loop splitting amplitudes for the dou- 
ble coUinear limit a ^ oi -|- 02 are known |5IHI7| . 
As a first step beyond the double coUinear limit, 
we have considered ^J the triple coUinear split- 
ting process q -^ q + q' + q', where q and q' de- 
note quarks of different flavours. To evaluate the 
one-loop splitting matrix we have used a process- 
independent method |4I7I1H| . Considering physi- 
cal spin polarizations, the splitting matrix is cal- 
culated from the sole Feynman diagrams where 
the parent parton emits and absorbs coUinear ra- 
diation. In the case q -^ q + q' + q' , one has to 
consider, for instance, the one-loop diagram de- 
picted in Fig.Efa). 

Our computation of Sp^ ' requires the evalua- 
tion of a set of basic one-loop (scalar and tensor) 
integrals. Besides the customary one- loop inte- 
grals, new integrals with additional propagators 
of the type l/{n ■ q) {q is the loop momentum and 
n is an auxiliary light-like vector), which come 
from the physical polarizations of the virtual glu- 
ons, have to be calculated. Some of these inte- 
grals, which resemble those encountered in axial- 
gauge calculations, were evaluated in Ref. in 
the context of the calculation of the one-loop 
double coUinear splitting a — > oi -f 02. More 
complicated integrals (higher-point functions) of 
this type are involved in triple coUinear split- 
ting processes. We have computed (to high or- 
ders in the e expansion) all the basic one-loop 
integrals that appear in any triple coUinear split- 
ting. These results can be applied to evaluate the 
one-loop splitting matrix of any splitting process 

a — > fli -I- 02 -I- 03 jni- 

The explicit expressions up to C(e*') of the 
splitting matrix Sp^^iq, and of the corresponding 

splitting function {PqJ„,) (see Fig.EIa)) are pre- 
sented in Ref. ^J . It is important to observe that 
{PqJqi ) has a contribution (which is proportional 
to the color factor (P^'^dabc) that changes sign by 
exchanging the momenta of the evolved quark and 



antiquark q' and q' . This charge asymmetry is a 
new quantum effect produced by the exchange of 
three gluons in the i-channel. When the charge 
asymmetry of {P^J„,) is combined with the cor- 
responding tree-level contribution (see Fig.|2Ib)), 
it leads to a non- vanishing value of the NNLO AP 
kernel Pqq/ — Pqqi . The main physical consequence 
of this effect is discussed in Sect. 0] 

4. STRANGE-QUARK ASYMMETRY 
IN THE NUCLEON 

Strange quarks and antiquarks play a funda- 
mental role in the structure of the nucleon |2]. 
Among the various strangeness-related proper- 
ties of the nucleon, the strange "asymmetry" , 
s{x) — s{x), in the densities of strange quarks 
and antiquarks, x being the light-cone momen- 
tum fraction they carry, is of particular interest. 
Since the nucleon does not carry any strangeness 
quantum number, the integral of the asymmetry 
over all values of x has to vanish: 



dx [s{x) — s{x)] — . 



(6) 



However, there is no symmetry that would pre- 
vent the X dependences of the functions s{x) and 
s{x) from being different. Therefore one can ex- 
pect s{x) y^ s{x), in general. 

Strange-antistrange asymmetries have been ex- 
tensively discussed in the literature. Various non- 
perturbative models of the nucleon structure |22| 
predict a fairly small value of the second moment 
of the strange-antistrange distribution, \{x{s — 
s))\ ^ 10"**. A global analysis of unpolarized 
parton distributions "2'6l reported improvements 
in the data fit if the asymmetry s{x) — s{x) is pos- 
itive at high X. However, a recent update |23] of 
this analysis reduces the asymmetry significantly. 
The most recent global QCD fit |2S| finds a large 
uncertainty for that asymmetry and quotes a 
range -0.001 < {x{s - s)) < 0.004. The strange 
asymmetry in the nucleon has become particu- 
larly relevant in view of the "anomaly" seen by 
the NuTeV collaboration in their measurement of 
the Weinberg angle 26^. The anomaly could be 
partly explained |27I28| by a positive value of the 
second moment {x{s — s)). 

The discussion reported so far regards strange- 
antistrange asymmetries that are generated by 
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Figure 3. Example of (a) virtual and (b) real 

(2)S 

contributions to Pns ■ 

non-perturbative mechanisms. Then, because of 
the customary scaling violation, the asymmetry 
becomes dependent on the hard-scattering scale 
Q at which the nucleon is probed. 

Perturbative QCD alone, however, definitely 
predicts a non-vanishing and Q-dependent value 
of the strange-antistrange asymmetry ^H]- The 
effect arises because at NNLO in perturbation 
theory the probability of the inclusive coUinear 
splitting (evolution) q -^ q' becomes different 
from that of q — > q' , and because the nucleon 
has u and d valence densities. 

Owing to charge conjugation invariance and 
flavour symmetry of QCD, the AP kernels that 
control the parton evolution of quark-antiquark 
asymmetries can be written as (see e.g. Ref. |29| ') 



9i9j 



p - — 



P- - — P- 



= S,, P^-^ + {Pqq, ~ Pqq,) . (7) 



The AP kernels Pqq, and Pqq, describe splittings 
in which the flavor of the quark changes, and 
Pqq' + Pqq' Starting from NNLO 12323 ■ In 
particular, as discussed at the end of Sect. O 
Pqq' -Pqq' = [as / (^tt))^ P^J '^ / N f , bccause of ths 
charge asymmetry produced by quantum effects 
at order ctg. The explicit expression of Pns is 
now available thanks to the recent computation 
by Moch, Vermaseren and Vogt 16 . 

The solution of the AP equations for the evo- 
lution (between the scales Qo and Q) of the A^- 



moment, (s — s) 



N 



„N-l 



(s — s)), of the strange- 



quark asymmetry reads ^Sl 

(s - s)^ (Q2) ^ Un{Q, Qo) [ (s - s)^^ (Ql) 



+ SP'^''{al{Q)~al{Qo))q'^\Ql) 



(V), 



(8) 



where q^^' = X]i=i(9i ~ ft) i^ the valence den- 
sity of the nucleon, Un is the evolution operator 
controlled by p(^\ and SPj^ is proportional to 



the NNLO kernel P^V^ HH]. At LO and NLO, 

(2) 

(5PjY = 0, and thus any asymmetry can only be 
produced by a corresponding asymmetry at the 
scale Qq. Starting from NNLO, the degeneracy of 
Pqqi and Pqqi IS removcd, and perturbative QCD 
necessarily predicts a non-vanishing s — s asym- 
metry driven by the valence density. 

Predictions for (s — s){x, Q^) based on Eq. ||S} 
have been presented in Ref. |18| . The asymme- 
try, s — s, is set to zero at a given low scale Qq, 
and then evolved upwards. The generated asym- 
metry is fairly sizable and turns out to be pos- 
itive at small x and negative at large x. Using 
Qo ~ 0.5 GeV (as in the 'radiative' parton model 
analysis of Ref. [33), a negative second moment 
is found: 

{x{s - s)) « -5 X lO-"* (g^ = 20 GeV^) . (9) 
The analysis has also be extended JSI to predict 
the asymmetries of heavy flavours c and b. 

Acknowledgments: G.R. thanks S. Moch and 
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